Polaron, molecule and pairing in one-dimensional spin-1/2 Fermi gas with an 

attractive Delta-function interaction 



O 
0-) 



C/2 



o 



ON' 

o; 



X: 



Xiwen Guan 

Department of Theoretical Physics, Research School of Physics and Engineering, 
Australian National University, Canberra ACT 0200, Australia 
(Dated: February 15, 2012) 

Using solutions of the discrete Bathe ansatz equations, we study in detail the quantum impu- 
rity problem of a spin-down fermion immersed into a fully ploarized spin-up Fermi sea with weak 
attraction. We prove that this impurity fermion in the one-dimensional (ID) fermionic medium 
behaves like a polaron for weak attraction. However, as the attraction grows, the spin-down fermion 
binds with one spin-up fermion from the fully-polarized medium to form a tightly bond molecule. 
Thus it is seen that the system undergos a cross-over from a mean field polaron-like nature into a 
mixture of excess fermions and a bosonic molecule as the attraction changes from weak attraction 
into strong attraction. This polaron-molecule cross-over is universal in ID many-body systems of 
interacting fermions. In thermodynamic limit, we further study the relationship between the Fred- 
holm equations for the ID spin-1/2 Fermi gas with weakly repulsive and attractive delta-function 
interactions. 
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I. INTRODUCTION 

The study of one-dimensional (ID) spin-1/2 Delta- 
function interacting Fermi gas [l|, [2| is an active area 
of research in the field of cold atoms. The fundamen- 
tal physics of the model with arbitrary spin population 
imbalance are determined by the set_of transcendental 
equations which were found by Yang l] using the Bethe 
ansatz hypothesis in 1967. The model displays a re- 
markable Fulde-Ferrell-Larkin-Ovchinnikov (FFLO)-like 
pairing [3-^, quantum phase transitions and quantum 
critical phenomena |9l-[ll|. It has a novel phase dia- 
gram caused by a difference in the number of spin-up 
and spin-down atoms (T2l - [l6| . The key feature of the 
phase diagram was experimentally confirmed by Liao et 
al [Ol at Rice University in the strongly attractive regime 
of fermionic ^Li atoms confined to the two lowest sub- 
hyperfine states. 

In 3D, for weak attraction limit, a spin-down fermion 
propagates almost freely in a spin-up medium. As at- 
traction increases the spin-down atom is dressed with the 
localized cloud of scattered surrounding fermions consti- 
tuting the Fermi polaron [T8l - [20| . However, if one consid- 
ers a small portion of spin-down fermions immersed into 
a fully polarized spin-up medium with strong attraction, 
it is seen that the system undergos a phase transition 
from a fully-polarized Fermi gas into a mixture of excess 
fermions and bosonic molecules [U, [l^] . 

Using variational ansatz Parish showed that there does 
not exist a true polaron-molecule transition in the ID 
highly polarized Fermionic system (23j . In fact, this vari- 
ational ansatz is not valid to capture the quasiparticle be- 
haviour for the ID interacting Fermi gas because it gives 
a divergent integral. The polaron-like effect can persist 
in the ID many-body system (23 |. In the present paper, 
using asymptotic solution of the Bethe ansatz equations, 
we present an analytical study of the quantum impu- 



rity problem in ID Fermionic medium. We prove that a 
spin-down fermion immersed into fully polarized spin-up 
medium with weak attraction is dressed with surround- 
ing fermions and behaves like a Fermi polaron, also see 
[25| . The spin-down fermions receive a mean field from 
the fully-polarized Fermi sea. In this limit, decoupling 
the two spin components gives a polaron-like quasipar- 
ticle associating with the spin-down fermion dressed by 
the particle-hole excitations. The mean field binding en- 
ergy and the effective mass of the polaron can be ana- 
lytically calculated from the discrete Bethe ansatz equa- 
tions. However, as an attractive interaction grows, the 
spin-down fermion binds with one spin-up fermion from 
the medium to form a tightly bound molecule. The cross- 
over is evidenced by the changes from a mean field at- 
tractive binding energy of the polaron with an effective 
mass m* = m to the binding energy of the single molecule 
with an effective mass m* — 2m as the attraction grows 
from c = to c = —oo. Here m is the actual mass of 
the fermions, and c is the interaction strength. Thus the 
system undergos a cross-over from a mean field polaronic 
nature into a mixture of excess fermions and a bosonic 
molecule as the attraction changes from weak into strong 
attractions. Furthermore, we discuss the relationship be- 
tween the Fredholm equations for the ID spin-1/2 Fermi 
gas with weakly repulsive and attractive Delta-function 
interactions. 



II. THE MODEL 



The model Hamiltonian [l|, 
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describes ID 5-function interacting spin-i Fermi gas of N 
fermions with mass m by periodic boundary conditions to 
a line of length L. The field operators (f>i and (j)^ describe 
the fermionic atoms in the states | f) and |) , respectively. 
We use units of h — 2m = 1 and denote coupling con- 
stant giD — lT?c/m with c = — 2/aiD where aio is the 

effective ID scattering length [llj aiD — + Aa±. 

Here asD is the 3D scattering length, a± ~ ^/h/(rruJ±) 
is the transverse oscillator length, and A « 1.0326 is a 
numerical constant. For repulsive fermions, c > and 
for attractive fermions, c < 0. 



For an irreducible representation R 



[27j . where N-f and N^^ are the numbers of fermions at 
the two hyperfine levels | t) and | l) such that Nf > 
N^^. The energy eigenspectrum is given in terms of the 

quasimomenta {ki} of the fermions via E = ^ J2jLi : 
satisfying the BA equations 



M- 



exp(ifcii) ]^ 



fcj - Aq - ic' ' 



N 

n 



Ml 

n 



i=i 
i = 1,2, 



,,7V, 



a = l,2,. 



A^ - ic 
.,Mi 



(2) 
(3) 



Giraud and Combescot [sTJ in the context of Fermi po- 
larons. The polaron-like effect in the ID Fermi-Hubbard 
model was studied by the variational anastz in . Here 
we prove that a single spin-down fermion immersed into 
the a fully ploarized spin-up Fermi sea with weak attrac- 
tion is likely to behave like a polaron in such fermionic 
medium. For the weak coupling limit L\c\ <^ 1, we find 
that either the spin-down fermion and a spin-up fermion 
from the medium form a pair fcj^^^ = p ± i/3 or the spin- 
down fermion with a quasimomentum — p propagates 
in the medium. Whereas the rest are N — 2 real roots 
{h} with i = 1, . . . ,iV - 2. 

We consider the quasimomenta of a pair k^^-f — p ±ij3 
and N ~ 2 real roots {fcj} with i — 1, . . . , — 2. From 
the discrete Bethe ansatz equations we have 
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with i — 1, . . . , N ~ 2. From the second equation in 
we determine the imaginary part j3 from 



with the quantum number Mi = and a notation 
c' = c/2. The parameters {Aq} are the rapidities for 
the internal hyperfine spin degrees of freedom. 

In thermodynamic limits, and for attractive regime, i.e. 
c < 0, quasimomenta {hi} of the fermions with different 
spins form two-body bound states, i.e., ka = ^ i^^, 
accompanied by the real spin parameter Aa [H, [2§| . Here 
a = 1 , . . . , Ml . The excess fermions have real quasimo- 
menta {kj} with i = 1, . . . , — 2Mi. From these root 
patterns, the BA equations ([2]) and (jS]) become 
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We see that for weak coupling limit ji — > ^J\c\/ L whereas 
for strong coupling limit (3 — > |c|/2. From the Bethe 
ansatz equations (jS]), we obtain 
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In the above equations a ^ j3. 



Using the B A equations (jH]) and ([5]) , we find that in weak 
coupling limit the roots satisfy the following polynomial 
equations 
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III. POLARON-MOLECULE CROSSOVER 
A. Polaron-like state 

McGuire studied exact eigenvalue problem of A^ — 1 
Fermions of the same spin and one fermion of the oppo- 
site spin in 1965 and 1966 [s^. He calculated the en- 
ergy shift caused by this extra spin-down fermion. The 
highly polarized Fermi system was studied recently by 



with i = 1, . . . , A^ — 2. According to the Fermi statistics, 
here rii = ±1, ±2, . . . , (A^ — 2)/2 and Up is an integer. 
The ground state configuration corresponding to = 0. 
In the weak coupling limit, we have f3 = y^\c\/L so that 
p ^ (3. Thus we see that the bound pair is not essential 
in the limit L\c\ ^ 1. The key feature of the model in 
this limit is that the spin-down fermion receives a mean 
field from the fully-polarized Fermi sea. We consider the 
case Up ^ 0, i.e. p ^ for excitations. From Eq. ©, 
we can calculate the energy of the system with a single 
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spin-down fermion 
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In the above equations, symmetrization of the N — 2 
quansimomenta in the medium was considered, i.e. the 
real roots associating with N — 2 spin-up fermions can 
be symmetrized by ki w —kj up to the order of c. This 
is mainly because the quansimomenta of the spin-up 
fermions just have an order of |c| deviation from the ones 
of the free spin-up fermions. For p = 0, the result (fTT|) 
coincides with the ground state energy of the Fermi gas 
with one spin-down fermion given by McGuire (soj (on 
page 125). If we consider the excitations for the weakly 
bound pair in the surrounding fully-polarized Fermi sea, 
i.e. p ^ 0, we can find an explicit relation of the en- 
ergy depending on the total momentum of the system. 
Defining total momentum of the system q, thus in weak 
coupling limit, we find a relation between p and q as 
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Substituting (|12p into (ITT|) . the last term in (ITT]) is can- 
celled out. Then we obtain an energy shift 

E{q,N,N^ = l)-E^iN^,0) « ep-b + ^ (13) 

2m* 

that behaves like a quasi-particle polaron. Where 
E^{N^, 0) — -^jz^N^T^^ is the kinetic energy of spin- 
up fermions. It is interesting to note that the attractive 
mean field binding energy 
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depends on number of spin-up fermions and interaction 
strength [l^. In the above equation (|13|) . the polaron- 
like state with an effective mass m* w m(l + 0{<?)) that 
is almost the same as the actual mass of the fermions in 
the limit L\c\ ^ 1. This is consistent with the result in 
[30I [3l| . The last term of the equation ([20]) indicates that 
an attractive interaction always enhances the effective 
mass of the polaron. The addressed "binding energy" 
can be rewritten as 
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fp-b = 2^F\1\ 



(15) 



where the dimensionless interaction strength is defined 
7 = c/n. The fermi energy is ep 



^\n^^\ This 

2m 3 



binding energy indicates a mean field effect. This is a 
ID analog of Fermi polaron-like state resulted in from 
the weak attraction between the impurity and the fully 
polarized Fermionic medium. This polaron-like state also 
exists for weakly repulsive interaction, where the spin- 
down fermion experiences a repulsive mean field energy 
shift. 

The mean field polaron-like state occurs for a few spin- 
down fermions immersed into a fully polarized Fermionic 
sea. In thermodynamic limit and in weak coupling 
regime, i.e., cL/N ~ 1, the ground state is the BCS- 
like pairing state with a pairing correlation length larger 
than the average interparticle spacing. The correlation 
function for the single particle Green's function decays 
exponentially, i.e., (■01 gV'i.s) — > e~^/^ with ^ = vp/IS. 
and s =ti 4-1 whereas the singlet pair correlation function 
decays as a power of distance, i.e., (V'^j! -j-V'^J 4,V'i,tV'i.i) ^ 
. Here A is the energy gap, and the critical exponents 
^ and are both greater than zero. However, once the 
external field exceeds the critical value, the Cooper pairs 
are destroyed. Thus both of these correlation functions 
decay as a power of distance and the pairs lose their dom- 
inance, where, molecule and excess fermions form the po- 
larized FFLO pairing-like phase. In this phase the spacial 
oscillations of pairing correlation are caused by an imbal- 
ance in the densities of spin-up and spin-down fermions, 
i.e., — nj,, which gives rise to a mismatch in Fermi sur- 
faces between both species of fermions. In ID, the pair 
and spin correlations with the spacial oscillation signa- 
ture are a consequence of the backscattering for bound 
pairs and unpaired fermions (ssj . In next section, we will 
see that the polaronic signature is significantly different 
from the molecule state where single spin-down fermion 
and a spin-up fermion from the medium with a strong 
attraction form a tight bound molecule of two-atom. 



B. The molecule state 

In order to catch the signature of molecule state, we 
consider a spin-down fermion immersed into fully polar- 
ized spin-up medium with strong attraction, i.e. L\c\ 3> 
1. We assume the bound pair k^,^ = p ± i/3 and N ~ 2 
real roots {h} with i = 1, . . . , iV - 2. From ^ and © 
with an odd number N^^ we find the real roots 
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with rij ~ ±1. ±3, . . . , ±(A^t ^ !)■ The pair excitations 
correspond to p 7^ 0. Then we obtain the energy 
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In the above equation, 



(18) 



is the kinetic energy of N-^ spin-up fermions with one 
spin-down fermion. If we consider the pair excitations 
with total momentum q, we find the relation between 
the p and the total momentum of the system q, i.e. 
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Substituting ([T^ into (|17p , then we obtain an energy 
shift 
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that behaves hke a molecule with a binding energy 

sIti ~ "2^y ^ sItI 



(21) 



in the strong attractive regime i|c| ^ 1. In the above 
equations = ^sx^-^^"""^ is the kinetic energy of 
spin-up fermions. However, the effective mass of the 
molecule 
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is almost twice the actual mass of the fermions. In the 
above equations, the variation of the chemical potential 
A/i = n^TT^. We see clearly that the system undergos a 
cross-over from a mean field polaron-like nature into a 
mixture of excess fermions and a bosonic molecule as the 
attraction changes from a weak attraction into a strong 
attraction. 



IV. GROUND STATE ENERGY 



real ki. Without losing generality, we assume Mi is odd 
and TV is even. Substituting this root patterns into the 
BA equations ([2]) and we find the following equations 
to determine the positive roots {Aq} and the positive real 
quasimomenta {kj} by 
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1,2,...,-^. In the equation a — 13 is excluded. By 
iteration, we obtain the ground state energy for weakly 
attractive regime 
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If we define linear density n = N/L and the density 
of down-spin fermions = Mi/L. In thermaldynamic 
limit, then the ground state energy per length ([25]) be- 
comes 



(26) 



that agrees with the result given in [32| . The the ground 
state energy (^5]) is also valid for weakly repulsive in- 
teraction, i.e., for c > 0. This gives a mean field effect 
for the ID delta-function interacting Fermi gas in weak 
coupling regime. 



A. Weak attraction 



In order to see physical signature of the ground state 
energy at vanishing interaction strength, we first focus on 
weakly attractive interaction in which two fermions with 
spin-up and spin-down states form a weakly bound pair. 
In this regime, the weak bound pair is not stable because 
the kinetic energy of the pair is larger than the binding 
energy. In the limit c — > 0~, the unpaired fermions sit 
on two outer wings in the quasimomentum space due 
to the Fermi statistics. In this weak coupling limit, 
i.e. L\c\ « 1, the imaginary part of the pseudomo- 
menta for a BCS pair is proportional to \fc[Tj. Thus the 
bound state has a small binding energy Cb = T?\c\lmL 
that is proportional to |c|. For arbitrary polarization, 
the system is described by M\ weakly bound pairs with 
fcP w Aq ± \\J cjL and iV — 2M\ unpaired fermions with 



B. Strong attraction 

For strong attraction, i.e. L\c\ 3> 1, (or say c ^ kp\ 
the discrete BA equations ^ and ([3]) with the root pat- 
terns fcP — Aaibiic for pairs and fcj for unpaired fermions 
can be linearized. For even iV^, we obtain the momenta 
of tight bound pairs and excess fermions [32] 

with integers rij = ±1, ±3, . . . , ±(iV — 2N^ — 1) and ria = 
±1, 2, . . . , ±(iV|^ — 1). In this scenario the bound states 
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behave like hard-core bosons due to Fermi statistics. The 
per length ground state energy of the model with strong 
attraction and arbitrary polarization is given by 



E 

T 
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where nj. = iV^/L and the binding energy e\ — 
the effective energy for unpaired fermions and pairs 



and 
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(30) 



(31) 



From the energies ([50)1 and ([5T|). we see that the bound 
pairs have tails and the interfere with each other. But, it 
is impossible to separate the intermolecular forces from 
the interference between molecules and single fermions. 



V. RELATIONSHIP BETWEEN THE TWO 
SETS OF THE FREDHOLM EQUATIONS 

The fundamental physics of the model are determined 
by the set of transcendental equations which can be 
transformed to the generalised Fredholm equations in the 
thermodynamic limit. This transformation was found by 
Yang and Yang in series of papers on the study of spin 
XXZ model in 1966, see an insightful article by Yang 
[3^ . For repulsive interaction, the Bethe ansatz quasi- 
momenta {ki} are real, but all {Aq} are real only for the 
ground state. There are complex roots of Aq called spin 
strings for excited states. In the thermodynamic limit, 
i.e., L,N ^ 00, N/L is finite, the above Bethe ansatz 
equations ^ and ^ can be written as the generalized 
Fredholm equations 
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where the integration boundaries Bi , B2 are determined 
by 

n:=N/L= / ri{k)dk, Mi/L = / r2{k')dk' . 
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In the above equations, we denote the function 
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with c > for repulsive regime and c < for attractive 
regime. The ground state energy per unit length is given 



by 



E ^ 



k^ri{k)dk. 
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The functions r„i{k) denote the Bethe ansatz root distri- 
butions in parameter spaces, i.e. ri(fc) stands for quasi- 
momenta distribution function, whereas r2{k) are the dis- 
tribution functions for rapidity parameter A in the BA 
equations ^ and ([3]). The ground state energy and full 
phase diagram can be obtained by solving analytically 
the Fredholm equations. 

For attractive regime, i.e. c < 0, the BA equations ^ 
and (O become 
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where 9{x) ~ 2arctana;, and Ij = ~{N — 2Mi — 
l)/2, -{N - 2Mi - 3)/2, . . . , (iV - 2Mi - l)/2 and Jj = 
-(Ml - l)/2, . . . , (Ml - 3)/2, (Ml - l)/2. In thermo- 
dynamic limit, we introduced the density of unpaired 
fermions pi{k) = dIj{k)/Ldk and the density of pairs 
P2{k) = dJj{k)/Ldk. They satisfy the following Fred- 
holm equations [H, |2^ 
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The linear densities are defined by 
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2 / P2{k)dk 
rQ2 

P2{k)dk. 



pi{k)dk, 



(41) 



The ground state energy per length is given by 



E : 



Q2 



-Q2 



{2k^ - c^/2) p2ik)dk+ / k^pi{k)dk. (42) 
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We will investigate the relationship between the two sets 
of the Fredholm equations for the Fermi gas with repul- 
sive and attractive delta-function interactions. 

In the light of Takahashi's unification of the ground 
state energy of the spin-1/2 weakly interacting Fermi gas 
[35I , we first examine the relationship between the Fred- 
holm equations for ID Fermions with repulsive and with 
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attractive delta-function interactions. It is convenient 
to use Yang's operator notations [11] for the Fredholm 
equations. Here we denote the integral operator 



Furthermore, in repulsive regime, the Fredholm equa- 
tions ([5^ and ([55)1 with the Fermi boundaries conditions 
(|Tf)) and exhibit a symmetry 
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Ai i — > A2, c 



Ai < — > A2, 
— > — c. 



(54) 



which is a symmetric function. We define the projection 
operators Ai and its dual projection operators Ai 



{k\A,n) = 0, for |fc| > B„ 

{k\A,n} =n{k), for |fc| < B„ 

{k\Airi) ^ ri{k), for |fc| > Bi, 

{k\A,n) ^ 0, for |fc| < B„ 



(44) 
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where i — 1,2. Similar notations are carried out for 
attractive interaction regime. 

The Fredholm equations (j3D) and (1551) for repulsive 
interaction regime can be rewritten in terms of these op- 
erators 



This symmetry relates to the spin-up and spin-down re- 
versal symmetry of the model. This transformation maps 
the eigenstates with down-spin atoms and up-spin 
atoms one-to-one onto the eigenstates with up-spin 
atoms and down-spin atoms for the gas. Similarly, 
in attractive regime, the Fredholm equations and 
(HUl) with the Fermi boundaries conditions (|5ip and 15/ 
preserve the spin reversal symmetry 
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where the integration boundaries Bi satisfy the following 
conditions 



As a consequence of the mapping (|53p . the ground 
state energy smoothly connects at vanishing interaction 
strength. To see this point, we need to unify the ground 
state energy for weakly repulsive and weakly attractive 
regimes. For weakly repulsive regime, the ground state 
energy is given by 
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for repulsive regime. Here we denoted 
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Substituting (j46l) into the above equations, we obtain the 
ground state energy per length 



For attractive regime the Fredholm equations are 
rewritten as 
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where Qi and Q2 are determined by 

\k\A2P2)G+iQi,k)dk, (51) 
{k\Aipi)G-{k,Q2)dk. (52) 



where 
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We prove that under a mapping 

ri(fc) < > pi{k), r2{k) 



P2{k), 



(53) 



gy{x) = 1 - G+{Bi,x). 

For attractive regime, the ground state energy (j42|) is 
rewritten as 



the Fredholm equations Eq. (gS]) with (gT]), (gll) for 
repulsive regime and the Fredholm equations (j50p with 
(ISTI) . (|5^ for attractive regime are identical. In the above 
equations c > for repulsive interaction regime and c < 
for attractive interaction regime are implied. 



E 
T 



-Q2 



2k^ 



-Qi 



c 

" "2 
1 

2^ 



{k\A2P2)dk 



{k\K^\k'){k'\A2P2)dk' 



(58) 



dk. 
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Substituting the Eq. ([50)1 into the above equation (l58l) . 
we obtain the ground state energy for weakly attractive 
regime 



Ql 



1 



E = ^ + — H{k,Qi)dk 



Q2 



{k\Ki\k'){k'\Aipi)dk' 



k'\>C 



(59) 

H{k,Qi)dk. 



We see that the ground state of the gas with a weakly 
repulsive interaction and a weakly attractive interaction 
can be unified through (|57| and (|59)) . The ground state 
energy Eqs. ([57)1 and Eqs. ([CTl) can be calculated in a 
straight forward way. In weak coupling regime, it cov- 
ers the result (I^Bl) . also see [s^l- The energy smoothly 
connects at vanishing interaction streng th ( but not ana- 
lytically connects at c = 0). Takahashi [3a] proved that 
the energy is infinitely differentiable at c = for a real 
value of c. But the two sets of the Fredholm equations 
turn to be divergent in the region c — > iO, see a discussion 

in (na). 



VI. CONCLUSION 

In conclusion, we have studied the polaron-molecule 
crossover in the ID spin-1/2 Fermi gas with an attractive 
delta-function interaction. We have found that a spin- 
down fermion immersed into a fully ploarized spin-up 
Fermi sea with weak attraction is dressed to form a Fermi 
polaron-like quasiparticle in the ID fermionic medium. 
The spin-down fermion receives a mean field attraction 
from the fully-polarized Fermi sea. However, as the at- 
traction grows, the spin-down fermion binds with one 



spin-up fermion from the fully-polarized medium to form 
a tightly bond molecule. We have presented the mean 
field binding energy and effective mass of the polaron in 
the weak attraction limit and also presented the binding 
energy and effective mass of the molecule in strong at- 
traction regime. The system undergos a cross-over from 
a mean field polaron-like nature into a mixture of ex- 
cess fermions and a bosonic molecule as the attraction 
changes from a weak attraction into a strong attraction. 
The asymptotic solutions of the discrete Bethe ansatz 
equations provides insight into understanding the mean 
field nature of the polaron-like state and the novel pair- 
ing in the interacting Fermi gas. For both weak and 
strong coupling regimes, we have obtained the ground 
state energy from the Bethe ansatz roots of bound pairs 
and excess fermions, where they have the interfere with 
each other. Furthermore, we have proved that the two 
sets of the Fredholm equations for the ID spin-1/2 Fermi 
gas with repulsive and attractive delta-function interac- 
tions are identical. The result we obtained for weak and 
strong attractions opens to experimentally study such 
mean field nature of polaron-like state and molecule sig- 
nature in ID trapped cold atoms. Current experiment is 
capable of catching the polaronic signature of ID inter- 
acting quantum gases of cold atoms. This can be possibly 
achieved by using a species selective dipole potential, see 
a recent experiment on the polaronic dynamics of the ID 
Bose gas |3q |. 
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